Adaptive introgression reveals the genetic basis of a sexually selected syndrome in wall lizards.

The joint expression of particular colors, morphologies, and behaviors is a common feature of adaptation, but the genetic basis for such "phenotypic syndromes" remains poorly understood. Here, we identified a complex genetic architecture associated with a sexually selected syndrome in common wall lizards, by capitalizing on the adaptive introgression of coloration and morphology into a distantly related lineage. Consistent with the hypothesis that the evolution of phenotypic syndromes in vertebrates is facilitated by developmental linkage through neural crest cells, most of the genes associated with the syndrome are involved in neural crest cell regulation. A major locus was a ~400-kb region, characterized by standing structural genetic variation and previously implied in the evolutionary innovation of coloration and beak size in birds. We conclude that features of the developmental and genetic architecture contribute to maintaining trait integration, facilitating the extensive and rapid introgressive spread of suites of sexually selected characters.

Twenty years on from Developmental Plasticity and Evolution: middle-range theories and how to test them.

In Developmental Plasticity and Evolution, Mary-Jane West-Eberhard argued that the developmental mechanisms that enable organisms to respond to their environment are fundamental causes of adaptation and diversification. Twenty years after publication of this book, this once so highly controversial claim appears to have been assimilated by a wealth of studies on 'plasticity-led' evolution. However, we suggest that the role of development in explanations for adaptive evolution remains underappreciated in this body of work. By combining concepts of evolvability from evolutionary developmental biology and quantitative genetics, we outline a framework that is more appropriate to identify developmental causes of adaptive evolution. This framework demonstrates how experimental and comparative developmental biology and physiology can be leveraged to put the role of plasticity in evolution to the test.

Maternal behavioral thermoregulation facilitated evolutionary transitions from egg laying to live birth.

Live birth is a key innovation that has evolved from egg-laying ancestors over 100 times in reptiles. However, egg-laying lizards and snakes can have preferred body temperatures that are lethal to developing embryos, which should select against prolonged egg retention. Here, we demonstrate that thermal mismatches between mothers and offspring are widespread across the squamate phylogeny. This mismatch is resolved by gravid females adjusting their body temperature towards the thermal optimum of their embryos. We find that the same response occurs in both live-bearing and egg-laying species, despite the latter only retaining embryos during the early stages of development. Importantly, phylogenetic reconstructions suggest this thermoregulatory behavior in gravid females evolved in egg-laying species prior to the evolution of live birth. Maternal thermoregulatory behavior, therefore, bypasses the constraints imposed by a slowly evolving thermal physiology and has likely been a key facilitator in the repeated transition to live birth.

The relative impact of parental and current environment on plant transcriptomes depends on type of stress and genotype.

Through developmental plasticity, an individual organism integrates influences from its immediate environment with those due to the environment of its parents. While both effects on phenotypes are well documented, their relative impact has been little studied in natural systems, especially at the level of gene expression. We examined this issue in four genotypes of the annual plant Persicaria maculosa by varying two key resources-light and soil moisture-in both generations. Transcriptomic analyses showed that the relative effects of parent and offspring environment on gene expression (i.e. the number of differentially expressed transcripts, DETs) varied both for the two types of resource stress and among genotypes. For light, immediate environment induced more DETs than parental environment for all genotypes, although the precise proportion of parental versus immediate DETs varied among genotypes. By contrast, the relative effect of soil moisture varied dramatically among genotypes, from 8-fold more DETs due to parental than immediate conditions to 10-fold fewer. These findings provide evidence at the transcriptomic level that the relative impacts of parental and immediate environment on the developing organism may depend on the environmental factor and vary strongly among genotypes, providing potential for the interplay of these developmental influences to evolve.

Chromosome-level genome assembly of Lilford's wall lizard, Podarcis lilfordi (Günther, 1874) from the Balearic Islands (Spain).

The Mediterranean lizard Podarcis lilfordi is an emblematic species of the Balearic Islands. The extensive phenotypic diversity among extant isolated populations makes the species a great insular model system for eco-evolutionary studies, as well as a challenging target for conservation management plans. Here we report the first high-quality chromosome-level assembly and annotation of the P. lilfordi genome, along with its mitogenome, based on a mixed sequencing strategy (10X Genomics linked reads, Oxford Nanopore Technologies long reads and Hi-C scaffolding) coupled with extensive transcriptomic data (Illumina and PacBio). The genome assembly (1.5 Gb) is highly contiguous (N50 = 90 Mb) and complete, with 99% of the sequence assigned to candidate chromosomal sequences and >97% gene completeness. We annotated a total of 25,663 protein-coding genes translating into 38,615 proteins. Comparison to the genome of the related species Podarcis muralis revealed substantial similarity in genome size, annotation metrics, repeat content, and a strong collinearity, despite their evolutionary distance (~18-20 MYA). This genome expands the repertoire of available reptilian genomes and will facilitate the exploration of the molecular and evolutionary processes underlying the extraordinary phenotypic diversity of this insular species, while providing a critical resource for conservation genomics.

Chromosome-scale genome assembly of the brown anole (Anolis sagrei), an emerging model species.

Rapid technological improvements are democratizing access to high quality, chromosome-scale genome assemblies. No longer the domain of only the most highly studied model organisms, now non-traditional and emerging model species can be genome-enabled using a combination of sequencing technologies and assembly software. Consequently, old ideas built on sparse sampling across the tree of life have recently been amended in the face of genomic data drawn from a growing number of high-quality reference genomes. Arguably the most valuable are those long-studied species for which much is already known about their biology; what many term emerging model species. Here, we report a highly complete chromosome-scale genome assembly for the brown anole, Anolis sagrei - a lizard species widely studied across a variety of disciplines and for which a high-quality reference genome was long overdue. This assembly exceeds the vast majority of existing reptile and snake genomes in contiguity (N50 = 253.6 Mb) and annotation completeness. Through the analysis of this genome and population resequence data, we examine the history of repetitive element accumulation, identify the X chromosome, and propose a hypothesis for the evolutionary history of fusions between autosomes and the X that led to the sex chromosomes of A. sagrei.

Environmentally induced DNA methylation is inherited across generations in an aquatic keystone species.

Transgenerational inheritance of environmentally induced epigenetic marks can have significant impacts on eco-evolutionary dynamics, but the phenomenon remains controversial in ecological model systems. We used whole-genome bisulfite sequencing of individual water fleas (Daphnia magna) to assess whether environmentally induced DNA methylation is transgenerationally inherited. Genetically identical females were exposed to one of three natural stressors, or a de-methylating drug, and their offspring were propagated clonally for four generations under control conditions. We identified between 70 and 225 differentially methylated CpG positions (DMPs) in F1 individuals whose mothers were exposed to a natural stressor. Roughly half of these environmentally induced DMPs persisted until generation F4. In contrast, treatment with the drug demonstrated that pervasive hypomethylation upon exposure is reset almost completely after one generation. These results suggest that environmentally induced DNA methylation is non-random and stably inherited across generations in Daphnia, making epigenetic inheritance a putative factor in the eco-evolutionary dynamics of freshwater communities.

A single locus regulates a female-limited color pattern polymorphism in a reptile.

Animal coloration is often expressed in periodic patterns that can arise from differential cell migration, yet how these processes are regulated remains elusive. We show that a female-limited polymorphism in dorsal patterning (diamond/chevron) in the brown anole is controlled by a single Mendelian locus. This locus contains the gene CCDC170 that is adjacent to, and coexpressed with, the Estrogen receptor-1 gene, explaining why the polymorphism is female limited. CCDC170 is an organizer of the Golgi-microtubule network underlying a cell's ability to migrate, and the two segregating alleles encode structurally different proteins. Our agent-based modeling of skin development demonstrates that, in principle, a change in cell migratory behaviors is sufficient to switch between the two morphs. These results suggest that CCDC170 might have been co-opted as a switch between color patterning morphs, likely by modulating cell migratory behaviors.

Unravelling the origin of the common wall lizards (Podarcismuralis) in south-eastern Europe using mitochondrial evidence.

The origin of the common wall lizards (Podarcismuralis) populations in south-eastern Europe (namely in Bulgaria and Romania), representing the north-eastern range border of this species, was addressed using mitochondrial DNA. We compared cytochrome b sequences from Bulgaria and Romania with those from the contiguous range in Central Europe that are available from previous studies. We recorded five main haplogroups in Bulgaria and Romania, belonging to the Central Balkan clade. However, haplogroup III was recorded in more localities than previously found. Additionally, signs of haplotype admixture were identified in several populations along the Danube River. The presence of the Southern Alps haplotype in one population from Otopeni, Bucharest (Romania) and its close phylogenetic relationships to north Italy populations suggests human-mediated introductions of this wall lizard clade in Romania. Our results confirm that P.muralis can have non-native lineages and admixture through active human-mediated transport.

Population Genomics of Wall Lizards Reflects the Dynamic History of the Mediterranean Basin.

The Mediterranean Basin has experienced extensive change in geology and climate over the past six million years. Yet, the relative importance of key geological events for the distribution and genetic structure of the Mediterranean fauna remains poorly understood. Here, we use population genomic and phylogenomic analyses to establish the evolutionary history and genetic structure of common wall lizards (Podarcis muralis). This species is particularly informative because, in contrast to other Mediterranean lizards, it is widespread across the Iberian, Italian, and Balkan Peninsulas, and in extra-Mediterranean regions. We found strong support for six major lineages within P. muralis, which were largely discordant with the phylogenetic relationship of mitochondrial DNA. The most recent common ancestor of extant P. muralis was likely distributed in the Italian Peninsula, and experienced an "Out-of-Italy" expansion following the Messinian salinity crisis (∼5 Mya), resulting in the differentiation into the extant lineages on the Iberian, Italian, and Balkan Peninsulas. Introgression analysis revealed that both inter- and intraspecific gene flows have been pervasive throughout the evolutionary history of P. muralis. For example, the Southern Italy lineage has a hybrid origin, formed through admixture between the Central Italy lineage and an ancient lineage that was the sister to all other P. muralis. More recent genetic differentiation is associated with the onset of the Quaternary glaciations, which influenced population dynamics and genetic diversity of contemporary lineages. These results demonstrate the pervasive role of Mediterranean geology and climate for the evolutionary history and population genetic structure of extant species.

A highly conserved ontogenetic limb allometry and its evolutionary significance in the adaptive radiation of Anolis lizards.

Diversifications often proceed along highly conserved, evolutionary trajectories. These patterns of covariation arise in ontogeny, which raises the possibility that adaptive morphologies are biased towards trait covariations that resemble growth trajectories. Here, we test this prediction in the diverse clade of Anolis lizards by investigating the covariation of embryonic growth of 13 fore- and hindlimb bones in 15 species, and compare these to the evolutionary covariation of these limb bones across 267 Anolis species. Our results demonstrate that species differences in relative limb length are established already at hatching, and are resulting from both differential growth and differential sizes of cartilaginous anlagen. Multivariate analysis revealed that Antillean Anolis share a common ontogenetic allometry that is characterized by positive allometric growth of the long bones relative to metapodial and phalangeal bones. This major axis of ontogenetic allometry in limb bones deviated from the major axis of evolutionary allometry of the Antillean Anolis and the two clades of mainland Anolis lizards. These results demonstrate that the remarkable diversification of locomotor specialists in Anolis lizards are accessible through changes that are largely independent from ontogenetic growth trajectories, and therefore likely to be the result of modifications that manifest at the earliest stages of limb development.

Extensive introgression and mosaic genomes of Mediterranean endemic lizards.

The Mediterranean basin is a hotspot of biodiversity, fuelled by climatic oscillation and geological change over the past 20 million years. Wall lizards of the genus Podarcis are among the most abundant, diverse, and conspicuous Mediterranean fauna. Here, we unravel the remarkably entangled evolutionary history of wall lizards by sequencing genomes of 34 major lineages covering 26 species. We demonstrate an early (>11 MYA) separation into two clades centred on the Iberian and Balkan Peninsulas, and two clades of Mediterranean island endemics. Diversification within these clades was pronounced between 6.5-4.0 MYA, a period spanning the Messinian Salinity Crisis, during which the Mediterranean Sea nearly dried up before rapidly refilling. However, genetic exchange between lineages has been a pervasive feature throughout the entire history of wall lizards. This has resulted in a highly reticulated pattern of evolution across the group, characterised by mosaic genomes with major contributions from two or more parental taxa. These hybrid lineages gave rise to several of the extant species that are endemic to Mediterranean islands. The mosaic genomes of island endemics may have promoted their extraordinary adaptability and striking diversity in body size, shape and colouration, which have puzzled biologists for centuries.

Evolution of the locomotor skeleton in Anolis lizards reflects the interplay between ecological opportunity and phylogenetic inertia.

Anolis lizards originated in continental America but have colonized the Greater Antillean islands and recolonized the mainland, resulting in three major groups (Primary and Secondary Mainland and Greater Antillean). The adaptive radiation in the Greater Antilles has famously resulted in the repeated evolution of ecomorphs. Yet, it remains poorly understood to what extent this island radiation differs from diversification on the mainland. Here, we demonstrate that the evolutionary modularity between girdles and limbs is fundamentally different in the Greater Antillean and Primary Mainland Anolis. This is consistent with ecological opportunities on islands driving the adaptive radiation along distinct evolutionary trajectories. However, Greater Antillean Anolis share evolutionary modularity with the group that recolonized the mainland, demonstrating a persistent phylogenetic inertia. A comparison of these two groups support an increased morphological diversity and faster and more variable evolutionary rates on islands. These macroevolutionary trends of the locomotor skeleton in Anolis illustrate that ecological opportunities on islands can have lasting effects on morphological diversification.

Evolvability and evolutionary rescue.

The survival prospects of threatened species or populations can sometimes be improved by adaptive change. Such evolutionary rescue is particularly relevant when the threat comes from changing environments, or when long-term population persistence requires range expansion into new habitats. Conservation biologists are therefore often interested in whether or not populations or lineages show a disposition for adaptive evolution, that is, if they are evolvable. Here, we discuss four alternative perspectives that target different causes of evolvability and outline some of the key challenges those perspectives are designed to address. Standing genetic variation provides one familiar estimate of evolvability. Yet, the mere presence of genetic variation is often insufficient to predict if a population will adapt, or how it will adapt. The reason is that adaptive change not only depends on genetic variation, but also on the extent to which this genetic variation can be realized as adaptive phenotypic variation. This requires attention to developmental systems and how plasticity influences evolutionary potential. Finally, we discuss how a better understanding of the different factors that contribute to evolvability can be exploited in conservation practice.

Plasticity and evolutionary convergence in the locomotor skeleton of Greater Antillean Anolis lizards.

Plasticity can put evolution on repeat if development causes species to generate similar morphologies in similar environments. Anolis lizards offer the opportunity to put this role of developmental plasticity to the test. Following colonization of the four Greater Antillean islands, Anolis lizards independently and repeatedly evolved six ecomorphs adapted to manoeuvring different microhabitats. By quantifying the morphology of the locomotor skeleton of 95 species, we demonstrate that ecomorphs on different islands have diverged along similar trajectories. However, microhabitat-induced morphological plasticity differed between species and did not consistently improve individual locomotor performance. Consistent with this decoupling between morphological plasticity and locomotor performance, highly plastic features did not show greater evolvability, and plastic responses to microhabitat were poorly aligned with evolutionary divergence between ecomorphs. The locomotor skeleton of Anolis may have evolved within a subset of possible morphologies that are highly accessible through genetic change, enabling adaptive convergence independently of plasticity.

Spatial variation in gene flow across a hybrid zone reveals causes of reproductive isolation and asymmetric introgression in wall lizards.

Hybrid zones provide insights into the evolution of reproductive isolation. Sexual selection can contribute to the evolution of reproductive barriers, but it remains poorly understood how sexual traits impact gene flow in secondary contact. Here, we show that a recently evolved suite of sexual traits that function in male-male competition mediates gene flow between two lineages of wall lizards (Podarcis muralis). Gene flow was relatively low and asymmetric in the presence of exaggerated male morphology and coloration compared to when the lineages share the ancestral phenotype. Putative barrier loci were enriched in genomic regions that were highly differentiated between the two lineages and showed low concordance between the transects. The exception was a consistently low genetic exchange around ATXN1, a gene that modulates social behavior. We suggest that this gene may contribute to the male mate preferences that are known to cause lineage-assortative mating in this species. Although female choice modulates the degree of reproductive isolation in a variety of taxa, wall lizards demonstrate that both male-male competition and male mate choice can contribute to the extent of gene flow between lineages.

Vertical Transmission of a Nematode from Female Lizards to the Brains of Their Offspring.

Parasites have evolved a diversity of lifestyles that exploit the biology of their hosts. Some nematodes that parasitize mammals pass via the placenta or milk from one host to another. Similar cases of vertical transmission have never been reported in avian and nonavian reptiles, suggesting that egg laying may constrain the means of parasite transmission. However, here we report the first incidence of transovarial transmission of a previously undescribed nematode in an egg-laying amniote, the common wall lizard (Podarcis muralis). Nematodes enter the developing brain from the female ovary early in embryonic development. Infected lizard embryos develop normally and hatch with nematodes residing in their braincase. We present a morphological and molecular phylogenetic characterization of the nematode and suggest that particular features of lizard biology that are absent from birds and turtles facilitated the evolutionary origin of this novel life history.

Enhanced locomotor performance on familiar surfaces is uncoupled from morphological plasticity in Anolis lizards.

The radiation of Anolis lizards in the Caribbean is associated with a diversification of the functional match between morphology, habitat use, and locomotor performance. It has been hypothesized that the microhabitat a lizard is reared in can achieve a similar fit of form and function within a species. This predicts that plasticity in the locomotor apparatus is accompanied by changes in perching behavior or improved locomotor performance. To test this, we raised juveniles of two species (Anolis sagrei and Anolis carolinensis) on either broad or narrow surfaces and examined perching behavior and locomotor performance as well as the shape of the pectoral and pelvic girdles, limb length, and thickness of the long bones. Perching behavior was not affected by the habitat surface experienced during ontogeny. However, individuals raised on broad surfaces showed better locomotor performance on broad surfaces, and the magnitude of the effect was as large as the difference between the two species. Both species showed modifications of pectoral and pelvic shape, but only A. carolinensis developed longer limbs on broad surfaces. However, these morphological adjustments induced by physical activity did not explain why lizards raised on broad surfaces performed better. Thus, it appears that early-life experiences can affect both the morphology of the locomotor apparatus and locomotor performance in Anolis lizards, without the two being functionally connected.

Developmental plasticity and evolutionary explanations.

Developmental plasticity looks like a promising bridge between ecological and developmental perspectives on evolution. Yet, there is no consensus on whether plasticity is part of the explanation for adaptive evolution or an optional "add-on" to genes and natural selection. Here, we suggest that these differences in opinion are caused by differences in the simplifying assumptions, and particular idealizations, that enable evolutionary explanation. We outline why idealizations designed to explain evolution through natural selection prevent an understanding of the role of development, and vice versa. We show that representing plasticity as a reaction norm conforms with the idealizations of selective explanations, which can give the false impression that plasticity has no explanatory power for adaptive evolution. Finally, we use examples to illustrate why evolutionary explanations that include developmental plasticity may in fact be more satisfactory than explanations that solely refer to genes and natural selection.

Evolutionary lability in Hox cluster structure and gene expression in Anolis lizards.

Hox genes orchestrate development by patterning the embryonic axis. Vertebrate Hox genes are arranged in four compact clusters, and the spacing between genes is assumed to be crucial for their function. The genomes of squamate reptiles are unusually rich and variable in transposable elements (TEs), and it has been suggested that TE invasion is responsible for the Hox cluster expansion seen in snakes and lizards. Using de novo TE prediction on 17 genomes of lizards and snakes, I show that TE content of Hox clusters are generally 50% lower than genome-wide TE levels. However, two distantly related lizards of the species-rich genus Anolis have Hox clusters with a TE content that approaches genomic levels. The age distribution of TEs in Anolis lizards revealed that peaks of TE activity broadly coincide with speciation events. In accordance with theoretical models of Hox cluster regulation, I find that Anolis species with many TEs in their Hox clusters show aberrant Hox gene expression patterns, suggesting a functional link between TE accumulation and embryonic development. These results are consistent with the hypothesis that TEs play a role in developmental processes as well as in evolutionary diversifications.